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INTRODUCTION  AND  REVIEW  OF  LITERATURE 

Variation  in  responses  of  insects  to  a  stimulus  is 
unavoidable  as  every  entomologist  knows.   A  knowledge  of  the 
type  and  extent  of  this  variance  is  necessary  in  the  design 
and  analysis  of  experiments. 

Bliss  and  Beard  (1935)  in  a  study  of  variation  of 
response  to  gases  of  the  adult  milkweed  bug  (Oncopeltus 
fasciatus.  Dall. )  stated  that  there  can  be  two  sources  of 
variation  in  susceptibility.   The  first  he  called  "static" 
which  would  cause  certain  individuals  to  maintain  their 
susceptibility.   The  second,  "dynamic  variation",  may 
change  the  susceptibility  of  the  individual  very  quickly. 
These  two  types  of  variation  were  also  discussed  by  Clark 
(1933).   Dynamic  variations  may  be  caused  by  several  factors. 
Rai,  et  al.  (1956)  found  that  malathion  was  more  toxic  to 
houseflies  (Musca  domestica.  Linn.  )  as  the  temperature  was 
increased  from  63°  to  75°  P..   Davidow  and  Hagan  (1955) 
stated  that  the  toxicity  of  a  compound  may  vary  from  species 
to  species  of  insect  and  with  environmental  conditions. 
They  suggested  that  sex,  age,  weight,  and  state  of  nutrition 
may  affect  susceptibility.   Intraspecies  differences, 
differences  which  normally  occur  in  a  population,  may  have  a 
marked  effect  on  results  obtained  in  acute  toxicity  studies. 
Yamasaki  and  Ishli  (1957)  found  that  variation  in  suscepti- 
bility can  be  introduced  by  species  variation,  sex,  stage  of 
growth  as  well  as  temperature  and  moisture  conditions. 


Nagasawa  (1957)  studied  the  effects  of  diet  on  knockdown  and 

length  of  adult  life  of  houseflies.   He  also  found  that  male 

houseflles  were  more  susceptible  to  p,p'-  DDT  powder  than 

were  the  females.   Sun  and  Sun  (1953)  found  that  young  house 

flies  are  more  susceptible  than  older  one3.   Gunther,  et  al. , 

(1958),  in  studies  of  toxicities  of  l\l\.   compounds  related  to 

diphenylme thane,  found  variation  in  replications  which  were 

repeated  at  different  times.   "Agreement  among  replicates  in 

a  given  serios  of  tests  was  excellent,  but  variations  among 

the  averages  of  series  repeated  at  different  times  were 

considerable  and  amounted  to  as  much  as  600/*»  in  the  i»T^0 

values  and  1000$  in  the  DD^Q  values. "  The  confused  flour 

beetle,  (Tribolium  confusum  Div. ) ,  was  used  as  the  biological 

indicator. 

Finney  (1952)  described  different  types  of  bioassays. 

One  type,  "indirect  assay  based  upon  quantal  ('all  or  nothing') 

responses",  has  a  unique  characteristic.   The  variance  of  any 

per  cent  mortality  can  be  calculated  according  to  a  formula 

discussed  by  Wadley  (19i|.8).   The  variance  of  p=  £J*»  when  p 

1 
represents  the  proportion  affected,  q  equals  (l-p)  and  N  the 

number  of  individuals  in  the  test.   Wadley  (19^8)  stated 

that  this  variance  of  a  proportion  (a  )  can  be  increased  by 

P 

a  number  of  factors.   These  include  environmental  ciianges, 
but  the  most  important  are  real  changes  in  resistance  of 
successive  lots  of  insects.   Heal  differences  in  susceptibil- 
ity, over  and  above  tne  true  variance,  is  common  among  such 


groups. 

It  ia  apparent,  then,  that  one  source  of  variation,  the 

2 
variance  of  p  (a  )  is  proportional  to  p  (Finney,  1952). 

Hence,  lioulden,  (1956)  discussing  probit  analysis,  stated 

that  the  variance  of  Y,  the  dependent  variable,  at  all  levels 

of  X,  the  independent  variable,  is  not  equal.   He  stated  that 

2 

a     is  greatest  when  p=*q=o.5  and  decreases  to  a  minimum  when 

p=0.0,  p=1.0.   This  variance  must  be  considered  in  the 
analysis  of  experiments.   Pitting  a  regression  formula 
througn  data  involves  making  the  assumption  of  homoscedastic- 
ity  (Snedecor  and  Cochran,  1957).   Goulden  (1956)  stated 
that  this  is  definitely  not  the  case  in  quantal  type  assays. 
Several  methods  of  analysis,  disregarding  heterogeneity 
of  variance,  have  been  suggested.   Wadley  (I9I4.8)  stated  that 
the  analysis  of  variance  could  be  used  when  H   is  large  or 
when  means  of  several  percentages  are  used.   Beru.son  (1950) 
tabulated  LD^'s  calculated  by  eight  different  methods  and 
compared  them.   lie  stated  that  for  practical  importance 
differences  among  them  are  quite  small.   Finney  (1952) 
recognized  that  one  of  the  oldest  and  most  widely  used 
principles  for  the  estimation  of  unknown  parameters  is  tne 
method  of  least  squares.   In  this  metnod,  the  estimates  are 
taken  as  those  values  which  minimize  the  sum  of  squares  of 
the  differences  between  the  observation  and  the  corresponding 
theoretical  values  or  expectations....   For  quantal  responses, 
minimization  of  the  quantity  Z/~fc[ (p-P)£7  has  no  merits  as  a 


method  of  estimation:   it  takes  no  account  of  the  fact  that 
the  variance  of  p(o )  depends  on  p,  so  that  large  deviations 
from  p  can  occur  much  more  easily  in  the  region  of  p=q=0.5 
than  at  the  extremes.   If  the  method  is  modified  and  based 
upon  a  weighted  sum  of  squares,  with  the  reciprocals  of  the 
variances  as  weights,  the  quantity  to  be  minimized  becomes 

^  Z.  " — wZ J   »  *^e  lQas*  squares  method  is  tnen  identical 

with  that  of  minimum  chi  square  (a  ). 

Bartlett's  test  of  homogeneity,  which  tests  the 

2  2  2 

hypothesis:  a-^  ~   a     ~  * %  *    -  o     is  outlined  by  Snedecor  and 

Cochran  (1957)  for  both  equal  and  unequal  subclass  numbers. 
Another  such  test  is  in  general/^   -,=  Pfi  when  N  is  equal 
for  each  per  cent  response  (Wadley  19^8).   Snedecor  and 
Cochran  (1957)  gave  a  method  ofX   for  unequal  numbers.   If 
the  number  in  a  group  (N),  the  number  affected  (X),  the 
per  cent  affected  (p),  and  the  product  (pX)  are  tabulated, 

chi-square  is  100/~  "pP  ( iqqIp ) ~J   with  a-1  degrees  of 

freedom,  where  a  is  the  number  of  groups. 

Houseflies  are  known  to  vary  in  response  to  a  toxicant 
when  various  conditions  are  changed.   However,  men  the  age, 
sex,  food,  and  rearing  conditions  are  neld  relatively 
constant  otner  variables  may  cause  different  lots  of  house- 
flies  to  respond  differently.   Lots  represent  different 
groups  of  houseflies  which  are  of  the  same  genetic  population 
but  from  different  parents. 


Variation  from  lots  to  lots  was  tested  to  determine  if 
this  variance  was  greater  than  the  characteristic  variance 
of  quantal  responses.   Tne  knowledge  gained  may  be  used  in 
designing  and  evaluating  future  experiments  making  critical 
comparisons  between  insecticides. 

MATERIALS  AND  METHODS 

Materials 

K.UN  strain  housef  lies  (lviusca  domestics  L.  )  were  used 
in  this  study.   These  were  originally  obtained  from  Kansas 
University  and  have  been  maintained  in  various  laboratories 
without  exposure  to  any  insecticide  for  15  to  20  years. 

Plies  for  these  tests  were  reared  in  tne  laboratory 
according  to  CSMA  procedure,  except  for  slight  modifications, 
A  description  of  the  rearing  procedure  is  given  by  Rai,  et 
al.  (1956). 

Malathion  used  for  testing  was  obtained  from  the 
American  Cyanamid  Company,  Stamford,  Connecticut  and  was 
99.6  per  cent  pure.   Glass-distilled  ecetone  was  used  for 
dilutions  and  cleaning  glassware.   The  insecticide  (mala- 
thion and  acetone)  was  applied  by  a  micro-device  designed 
and  described  by  Roan  and  Maeda  (1953). 

The  houseflies  were  conditioned  and  held  during  the 
test  in  a  constant  te.nporature  cabinet.   The  temperature 
in  the  6.2  cubic  foot  box  was  held  at  76  t   2  degrees 
Fahrenheit  during  the  test.   Circulation  within  the  chamber 


was  provided  by  a  fan,   A  detailed  description  of  the  test 
procedure  was  described  by  Rai,  et  al.  (195>b). 

Methods 

A  complete  block  design  was  used  to  study  the  response 
of  male  houseflies  to  an  insecticide.   Days  were  used  as 
blocks  with  repetitions  of  six  levels  of  insecticide 
appearing  six  times  in  each  block.   The  levels  of  insecti- 
cide (expressed  in  ug./ml.  per  fly)  were  O.lj.0,  0.35#  0.30, 
0.25,  0.20,  and  0.00  with  the  0.00  level  as  a  measure  of 
natural  mortality.   Pour  day  old  male  houseflies  were  used 
in  the  tests. 

Shortly  before  treatment,  samples  of  malathion  were 
weighed  and  five  dilutions  prepared.   Test  solutions  were 
prepared  each  test  day. 

Each  day,  six  repetitions  of  £0  flies  each  were  tested 
at  each  level  of  insecticide.   A  control,  for  each  repeti- 
tion, which  received  no  insecticide  but  an  equal  exposure 
to  carbon  dioxide,  was  maintained.   In  this  manner  any 
mortality  due  to  factors  otner  than  the  insecticide  and 
acetone  could  be  tabulated. 

A  minimum  of  carbon  dioxide  was  used  to  anesthetize  the 
flies  while  sexing  and  treating.   All  parts  of  the  dilution 
and  treatment  apparatus  were  thoroughly  washed  and  rinsed 
with  acetone  before  and  after  their  use.   The  syringe  and 
needle  were  rinsed  with  glass-distilled  acetone  after 


completing  treatment  with  each  level.   They  were  then  dried 
with  a  stream  of  air.   The  syringe  was  rinsed  in  the  test 
solution,  then  filled  for  treatment.   The  first  few  drops 
of  insecticide  from  the  syringe  were  discarded  to  free  the 
syringe  from  air  bubbles. 

Individual  treatment  was  accomplished  by  placing  the 
needle  to  the  mesosternum  of  each  fly  and  applying  the 
insecticide.   The  treated  fly  was  immediately  placed  in  a 
recovery  jar  and  a  second  fly  treated.   When  50  flies  were 
in  a  jar,  cheesecloth  was  placed  over  the  mouth  of  the  jar, 
and  held  fast  by  a  rubber  band.   Recovery  jars  were  then 
placed  on  the  bottom  shelf  of  the  constant  temperature 
cabinet. 

Following  each  test  the  recovery  jars  were  rinsed  with 
acetone,  scrubbed  with  hot  detergent  water,  rinsed  with 
acetone  again  and  allowed  to  dry.   Before  reuse,  each  jar 
was  rinsed  with  acetone  a  third  time  and  allowed  to  dry. 
Cheesecloth  and  food  were  discarded  after  each  test. 

Plies  remained  in  the  constant  temperature  cabinet  for 
2k   hours  following  the  application  of  the  insecticide. 
They  were  then  removed  and  counts  of  surviving  insects 
made. 

The  average  weight  of  a  particular  group  of  houseflies 
was  obtained  by  weighing  a  random  group  of  j?0  flies.   This 
average  weight  was  used  to  calculate  the  LIVq  estimates  for 
each  repetition. 
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RESULTS 

This  experiment  was  begun  June  1,  1958.   It  was 
designed  to  measure  the  variance  of  a  bioassay  with  quantal 
response  when  repetitions  are  conducted  over  a  six  week 
period.   Actual  testing  was  done  every  weeit  from  June  5 
through  July  23. 

Response  was  measured  by  counting  fly  mortality  in  each 
jar  of  50  male  houseflies.   i^ach  day  six  repetitions  of  $0 
male  flies  at  six  levels  of  toxicant  were  performed.   One 
level  served  as  a  control. 


Table  1.   Average  per  cent  mortality  of  houseflies 
at  five  concentrations  of  malathion  and  a  control  on  six  dates. 


Day 

:Levels 

of 

mala 

thion 

in 
ace 

I  microf 
tone 

irams 

per  milliliter 

:'  .kO 

: 

• 
• 

• 

.35  : 

.30 

• 
• 

• 
• 

• 

.2S      1 

.20 

• 
• 

:  Control 

6-18-58 

83 

73 

k$ 

U5 

2$ 

.0 

6-25-58 

* 

95 

86 

83 

k6 

1.0 

7-  2-58 

95 

89 

83 

60 

16 

1.0 

7-  9-58 

98 

91 

78 

Til- 

& 

3.0 

7-16-58 

* 

89 

76 

53 

37 

.0 

7-23-58 

96 

91 

72 

59 

30 

.0 

Table  1  illustrates  that  there  are  differences  in 
mortalities  at  different  levels  of  toxicant.   This  is 
expected,  for  the  assigned  levels  were  pretested  to  insure 
covering  a  minimum  of  50  per  cent  of  the  mortality  range 


with  each  test.   It  is  also  evident  that  within  a  concentra- 
tion there  is  variation  in  response  from  day  to  day.   Control 
mortality  was  less  than  three  per  cent  on  all  days. 

It  was  found  in  the  literature  that  such  measures  us 
"dead  or  alive"  are  quantal  in  nature  and  have  a  variance 
which  is  subject  to  change  as  the  proportion  affected 
changes.   It  is  possible  to  test  for  differences  in  variances 
at  each  average  p  by  using  Bartlett's  test  of  homogeneity 
(Snedecor  and  Cochran,  1957).   Corrected  chi  square  by 
Bartlett's  test  was  3i|.  82. 

Chi  square  exceeds  the  tabled  value  at  P  =  .00$.   From 

2 
this  evidence  it  was  concluded  that  s  is  different  for 

different  p.   In  Table  3  is  given  the  average  per  cent 

mortality  at  each  level  of  toxicant.   HeterOf. eneity  of 

variance  would  be  expected  from  such  data,  for  the  variance 

of  p  (estimated  by  s^)  would  be  greatest  when  F  ■  0.$  and 

a  minimum  when  p  =  1  or  0.   This  is  obvious  when  p  ana  the 

2 

corresponding  s  are  compared  to  other  values. 

2 
It  is  also  noticed  that  s  is  generally  quite  large. 

2 

The  estimates  of  s  for  p  -  77,  62,  and  35  are  especially 

large. 

2 
It  is  apparent  that  the  estimates  of  s   are  biased 

upward  by  day  to  day  differences.   It  is  possible  to  separate 

the  day  to  day  effects  and  obtain  an  unbiased  estimate  of 

2 
s  ,  by  the  analysis  of  variance. 

A  one-way  analysis  of  variance  was  calculated  for  each 


level  of  insecticide.   By  a  variance  component  analysis 

(Table  2)  an  estimate  of  a     was  obtained,  a  similar  form 

e 

is  given  by  Snedecor  and  Cochran  (1957). 
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Table  2.   Component  analysis  for  one-way  analysis 
of  variance  with  degrees  of  freedom. 


Source  of  Variation  :  D/P 


Days 
Within  Days 


Expectation  of 
mean  square 


30      of 


Estimates 


2     A,2 

s   +  6k0 

e       R 


By  the  one-way  analysis  of  variance  each  average  p  has 
an  unbiased  estimate  of  s   (Table  3). 

P  P 

Table  3.   Comparison  of  s  ^and  sf  1<or  ;flve 
average  p  with  a   . 


p 

*a2 

P 

:              2 

:                     s 

• 
• 

i         -5 

• 

93 

17.0172 

kl.  2825 

21.5323 

88 

29. 0987 

73.61^5 

28. 8889 

77 

53.9589 

237.6730 

62. 5W+5 

62 

65.2570 

208.  2635 

ij.6. 0222 

35 

62. 8190 

213.6250 

59. 0222 

•     2 
a     = 
P 

i 

Scheffe  (1956),  and  Rama ch and ran  (1958)  gave  a 
2      (n-l)sd        /_/      o  _ 
n-1  =    a2     where  £  (n-l)s  _/  is  equal  to 
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X  (Xi-*): 

{•I 

2  2  2  2 

Using  a  as  an  estimate  of  a     and  35(s  )  as  2b-,  all  corres- 

P    r-  1 

ponding  X  '  values  are  greater  than  the  tabular  value  at 

P  =  .05.   This  is  evidence  at  P  =  .  05  that  s2  is  too  large 

to  be  called  an  estimate  of  a   .   There  is  additional 

P 

variability  included  in  the  estimate,  sL,  which  overshadows 

the  variance  of  p.  (a  ). 

P 

2  2 

If  a  similar  test,  using  s   as  an  estimate  of  s  ,  the 

e 


formula  becomes: 


i 


Scheffe  (1956)  stated  that  s  is  an  unbiased  estimate. 

6 

p 

Chi  square  values  calculated  on  the  basis  of  s   when 

e 

compared  to  tne  tabular  value  proved  to  be  non-significant 

at  P  =  .05.   There  is  evidence  that  the  effects  of  days 

p  p 

bias  the  estimate  of  s   upward,  since  at  P  =  ,05  Sg  is  a 

p 
reliable  estimate  of  a   ,  the  true  variance  of  p. 

P 
By  the  formula  (-3-),    the  variances  of  percentages 

differ  as  the  per  cent  (p)  takes  on  different  values. 

These  differences  are  shown  to  be  significant  by  Bartlett's 

test  of  homogeneity.   The  variance  of  the  per  cent  when 

p 

estimated  by  s  is  too  large  to  be  recognized  as  an 

2 
estimate  of  o  .   by  removing  the  variance  caused  by  day  to 

day  differences,  a  second  estimate  of  a     is  found,  s2. 

P  © 
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This  "within"  variance,  sf,  la  a  reliable  estimate  of  a^  as 

Jr 


e; 


shown  at  P  =  .05  with  the  chi-square  test.   It  was  apparent 
that  the  effects  of  day  to  day  differences  were  present  and 

of  important  magnitude. 

2  P 

It  has  been  shown  that  s0  is  an  estimate  of  a   . 

e  P 

Another  variance  component  remains,  k^,  which  is  the  variance 

which  estimates  the  day  to  day  variation  at  a  particular 

level  of  toxicant.   It  mi^ht  be  expected  that  tnis  estimate, 

2 
k,,  would  be  constant  at  each  level;  i.e.  the  effects  of 

days  could  be  measured  equally  well  at  any  level. 

Using  a  general  formula  given  by  Bross  (1950)  it  is 


shown  that  k 


Among,  Ms  -  Within  Ms 

\  :   where  \   is  a 

2 


coefficient  assigned  to  k,,  in  the  component  analysis.   Thus 

.2  _  Among  Ms  -  Vvithin  Ms    „  ,      ^  ,  2 

«jj  -  5 .   values  of  k^  are  shown  in 

Table  k» 


Table  If.   Estimates  of  k  for  five  averuge  per  cent 
mortalities  and  inxraclass  correlation. 


Average  per  cent  mortality 

93 

88 

77 
62 
35 


Intraclass  correlation 
I 


19.71 

.U779 

5k*$$ 

.6538 

20k.  22. 

.7656 

189.28 

.80i4f 

180. Ok 

.7531 

Following  the  metnod  given  by  .bross  (1950)  for 
setting  a  confidence  limit  on  k.f    it  was  found  that, 
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P(6,9532aSkf  -£83.56!|5)  -  .90.   Since  this  confidence 

interval  does  not  include  all  estimates  of  k  ,  it  is 

2  R 

apparent  the  k  is  not  a  constant  from  level  to  level.   For 

n 

some  reason  day  effects  are  greater  when  the  proportion 
affected  is  near  0.5.   This  is  beyond  the  scope  of  these 
investigations. 

Intraclass  correlations  calculated  for  each  average 
p,  show  that  good  reproducibility  exists  within  each  level 
for  each  lot  of  houseflies. 

Prom  Bartlett's  test  of  homogeneity  it  was  found  that 

the  veriance  of  p  changed,  as  expected,  at  different  levels 

of  toxicant.   Finney  (1952)  has  shown  that  a  method  of 

p 

analysis  which  minimizes  the  quantity   2  /""  ,n*P"F',   T 

u    pq    -' 

takes  account  of  different  variances  of  p.   Such  a  method 
of  analysis,  wnich  fits  a  straight  line  through  observed 
values  of  p  for  each  level  of  concentration,  is  termed 
prcbit  analysis. 

To  expedite  computation  the  probit  analysis  was  done 
by  the  IBM  6^0  Data  Processing  Machine.   The  program  used 
for  the  analysis  was  written  by  Dr.  R.  R.  Sokal  cf  Kansas 
University. 

Results  obtained  from  this  analysis  included  estimates 
of  Beta  and  the  LD^Q,  with  5  per  cent  fiducial  limits.   The 


1 

Received  by  personal  communication. 
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estimated  mortality  (Y. )  for  each  concentration  (X.)  was 
also  computed  so  that  graphing  could  easily  be  accomplished. 

By  coding  the  individual  doses  in  terms  of  micro-grams 
per  gram  body  weight,  tbe  computed  ED^0  values  are  actually 
LIV  values. 

The  probit  analysis  was  computed  for  each  individual 
repetition  and  the  mean  of  all  36  repetitions.   The  hypoth- 
esis:   =  0  was  tested  at  P  =  .0$,  as  was  heterogeneity  chi 
square. 

Of  special  interest  was  the  hD^p»  as  it  is  a  measure 
of  the  amount  of  toxicant  to  create  a  £0  per  cent  mortality 
among  the  individuals  tested.   It  is  a  measure  of  suscepti- 
bility or  the  ability  of  a  group  to  tolerate  the  toxicant. 

It  is  known  that  an  estimate  of  a  parameter  is  approxi- 
mately normally  distributed  (Mood,  1930)  and  such  measure- 
ments can  be  studied  by  the  analysis  of  variance.   By  using 
the  functional  analysis  of  variance  described  by  LeClerg(l957) 
it  can  be  determined  by  testing  orthogonal  polynomials  if 
a  graphic  relationship  can  be  used  to  explain  the  response. 

Table  5.   Functional  analysis  of  variance  of  DD[-n 
values  determined  on  six  equally  spaced  days.-3 

:     :         :         : Statistical 
Source  of  variation  :  D/F  :    MS    :   P {Significance 

Days  5  "23.90 

Linear  1  !+$'.  03 

Quadratic  1  25.76 

Cubic  1  19. 72 

Quartic  1  1+. 65 

Remainder  1  2k.55 

Within  Days 3_0 1.07 


22.36 

.01 

1+2.12 

.01 

2k.  10 

.01 

18.1|5 

.01 

4.35 

.05 

22.96 

.01 
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From  Table  5  it  can  be  seen  that  the  test  to  test  effects 
on  the  IiD^0  were  significant  at  P  =  .01.   This  leads  to  the 
rejection  of  the  H  of  no  difference,  and  the  acceptance  of 
the  H  that  the  day  to  day  effects  are  different.   An  attempt 
to  find  a  polynomial  which  would  describe  the  day  to  day 
effects  proved  futile.   After  fitting  the  quartic  polynomial, 
a  significant  remainder  was  left.   It  is  concluded  that  day 
to  day  variation  in  UVq  values  are  net  systematic.   How- 
ever some  tendency  for  a  graphic  trend  is  present  as  is 
indicated  by  significant  mean  squares  for  linear,  quadratic, 
and  cubic.   The  quartic  polynomial  accounted  for  very  little 
variation  (3.9$>)  even  though  its'  mean  square  was  signifi- 
cant at  ?  =  . 05. 

The  LIVq  estimate  of  June  18  was  the  greatest  and  is 
different  from  the  other  estimates.   The  mean  LD^n  for  all 
tests  is  21.19  ug/gram  fly.   P(19.  89  ^£  U)^-^  22.27)  =  .95 
Is  the  probability  statement  based  on  .05  fiducial  limits. 

Figure  1  illustrates  the  differences  among  extreme 
regression  for  six  repetitions  when  pooled.   The  LD^0  for 
June  25  is  11.63  ug/gram  fly,  while  2]+.  83  ug/gram  fly  was 
the  HVq  for  June  18.   They  are  compared  to  the  average  for 
all  tests. 

The  mean  LIW  estimate  for  all  tests,  21.19  ug/gram 
fly,  is  compared  to  other  pooled  estimates  in  Table  6. 
Variation  In  LDj-q  estimates  is  shown  by  comparing  them  to 
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21. 19  ug/gram  fly.   Estimates  ranged  from  $5  per  cent  to 
117  per  cent  of  the  mean  estimate. 

Table  6.   Pooled  LD^q  estimates  from  six  tests  compared 
to  the  LD?q  estimate  for  all  tests. 


Date  tests 
performed 

:      Pooled 

LD^q    es 

itimate 

:    Pooled  estimates   as 
:    a  per  cent   of  21.19 
:               uprV.rain  fly 

June   18 

2k.  83 

117.18 

25 

11.63 

5k.  88 

July     2 

21.08 

99.1f8 

9 

23.00 

108.5^ 

16 

22.15 

10J|.  53 

23 

16.14 

76.17 

Comparisons  among  the  estimate  of  21. 19  ug/gram  fly 
and  estimates  collected  by  other  workers  using  similar 
techniques  and  KUN  strain  males  are  shown  in  Table  7. 


Table 

7. 

Estimates 
malathion 

Of 

on 

LD<n  for 
housefllc 

pure 
i3    by 

(99. 6£) 

various  workers. 

Name 

: 
: 

• 
• 

Dates      ; 

IJVq  estimate 

Insect 

L.    Rai 

1956 

20.83 

Male  KUN  housefly 

G.    Ware 

1957 

16.89 

Male  KUN  housefly 

G.    Krause 

1958 

21. 19* 

Male  KUN  housefly 

contains  effects  of  acetone. 

Estimates  of  $    were  obtained  from  the  probit  analysis. 
Each  estimate  was  tested  at  P  a  ,05  by  a  t  test,  the  null 
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Pooled  probit  regression  of  tests  performed 
June  18,  1958  and  June  25,  1958  compared  to 
average  for  all  tests  (LD^q). 

.. _ , June  Z$,    1958. 

June  18,  1958. 

All  tests. 
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hypothesis  was  HQ:  ^  =  o.   In  all  cases  this  hypothesis  was 
rejected. 

All  tests  were  conducted  using  the  same  population  of 
houseflies.   To  determine  the  homogeneity  of  tolerances, 
the  Q  estimates  were  separated  by  the  analysis  of  variance 
and  Duncan's  new  multiple  range  test  (LeClerg,  1957). 
Graphic  trends  were  tested  by  using  the  functional  analysis 
of  variance  and  testing  individual  degrees  of  freedom. 


Table  8.   Analysis  of  variance  of  Beta  estimates 
with  individual  degrees  of  freedom. 


Source   of  variation 

D/F        : 

M.S. 

:        SIG 

Repetitions 

5 

11.3532 

.01 

Linear 

1 

4.7392 

ns 

Quadratic 

1 

3.2623 

ns 

Cubic 

1 

15.6271 

.01 

Quart I c 

1 

2.5981 

ns 

Other 

1 

30.1+823 

.01 

Yiithin  Repetitions 

30 

1.3302 

The  analysis  of  variance  shows  that  differences  exist 
among  estimates  of  Beta,  taken  on  different  days.   However, 
the  analysis  of  individual  degrees  of  freedom  shows  tnat 
time  trends  explain  only  a  small  per  cent  (1+6. 3#)  of  the 
variability  in  the  Beta  estimates. 
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Table  9.   Duncan* 3  new  multiple  range  test,  separation 

of  average    estimates. 

Day         6/18    7/9    6/25    7/l6    7/23    7/2 
Ranked  means  5.38    $\Z       6.17    6.51+    7.36    9.01 
Significance  ;__;>______ 


means  underlined  are  not  significantly  different  at  i>  =  .05. 

Duncan's  mean  separation  test  (Table  9)  shows  that  the 
estimate  of  July  2  was  higher  than  all  others.   This 
indicates  that  the  homogeneity  of  tolerance  among  the  test 
group  was  greatest  on  this  date.   Lower  Beta  estimates  on 
other  dates  indicate  heterogeneity  of  tolerance  within 
test  groups.   Estimates  collected  on  June  2S,   July  l6,  and  23 
were  not  significantly  different.   The  second  highest  esti- 
mate 7.36  was  different  from  the  two  lowest  estimates. 

DISCUSSION 

Variation  in  response  was  studied  by  several  statistical 
methods.   It  was  found  that  quantal  type  responses  have  a 
variance  which  is  characteristic  to  binomially  distributed 
data.   However,  other  variables  also  affect  the  response  of 
male  houseflies. 

Cnaracteristic  to  a  binomial  distribution  is  the  fact 
that  the  variance  of  p  is  greatest  when  it  is  near  0.5,  and 
diminishes  as  it  approaches  1  or  0.   Bartlett's  test  of 
homogeneity  shows  that  the  variance  of  different  percentages 
are  significantly  different  at  the  .005  level  of  probability. 
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A  chi-square  test  suggests  that  lot  to  lot  differences 

are  great  enough  to  prevent  the  recognition  of  the  true 

o 

variance  a   . 
P 

Differences  in  response  from  lot  to  lot  were  signifi- 
cant.  There  was  no  discernable  trend  of  response  over 
equally  spaced  days  which  would  describe  tne  cnange  in 
susceptibility.   UVq  estimates  from  pooled  data  from  six 
daily  repetitions,  compared  to  tlie  U>£q  estimate  for  36 
repetitions  range  from  54*9  to  117.18  per  cent  of  the  mean 
estimate. 

The  mean  LD^q  ©stimate  calculated  by  pooling  all  36 
repetitions  was  21.19  ug/gram  fly.   This  compares  with 
20.83  ug/gram  fly  found  by  Lallan  Rai  in  1956  using  similar 
techniques  and  KUN  strain  males. 

It  was  found  by  the  analysis  of  variance  and  a  component 
analysis  that  the  variance  due  to  days  (lots  of  houseflies) 
was  significant  at  P  =  .  0£.   Intraclass  correlations  show 
there  was  greater  reproducibility  within  days  than  among 
days. 

The  greatest  source  of  variation  beyond  that  of  the 
binomial,  was  that  of  k^.   The  binomial  variance  cannot  be 
minimized  beyond  ^,  therefore  sample  size  is  determined 
by  the  variance  of  days. 

Various  numbers  of  lots  to  be  sampled  to  be  within  a 
specified  difference  of  the  UVq  with  two  levels  of  proba- 


22. 


bility  are  given  in  Table  10.   It  is  noticed  that  N  becomes 
very  large  as  d  becomes  small. 

Table  10.   Probability  of  an  estimate  being  within  a 
per  cent  of  the  LD^q  (malathion)  when  sampling  various 
numbers  of  lots  of  KUN  male  houseflies. 


18.9 
16.5 

15.5 

14.3 

13.4 
3.0 


d  as  a  per  cent   • 

Lots  of  houseflies  (il) 

of  LD5Q 

P  »  .90       :     P  *  .95 

6 

I 

9 
10 

13 
18 

?° 

ko 

120 


8 

10 

n 

12 

30 

120 


% 


This  is  evidence  that  sampling  days  (li)    is  more  useful 
and  efficient  in  determining  a  reliable  estimate  of  the 
^50  than  lar8e  numbers  of  repetitions  within  days. 
Repetitions  within  days  should  be  limited  to  the  necessary 
number  required  to  obtain  an  estimate  of  the  response  on 
that  day. 

Relative  efficiency  was  calculated  according  to  a 
method  found  in  Cochran  and  Cox  (1957).   It  was  found  that 
a  randomized  complete  block  design,  with  days  as  blocks 
has  a  relative  efficiency  of  1. 96  compared  to  a  completely 
random  design.   This  may  be  interpreted  to  mean  that 
approximately  one  half  as  many  tests  need  to  be  run  to 
obtain  an  equally  reliable  estimate  when  the  effects  of 
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days  are  removed. 

Table  11.   Number  of  lots  of  houseflies  to  be  sampled  to 
be  within  d=lf  ug/gram  fly,  at  various  levels  of  oi   and  %  , 

P  of  type  I  error  :   P  of  type  II  error   :   Lots  to  detect 
©£ : 0 : d 

.05  .3  19 

.05  .2  25 

.01  .3  30 

.01  .2  36 

The  effect  of  power  of  the  test  (1-3)  upon  the  number 
of  lots  to  be  sampled  is  illustrated  in  Table  11.   Within 
each  level  of  0(  the  necessary  number  of  lots  to  detect  a 
difference  of  I4.  ug/gram  fly  is  greater  when  the  power  of 
the  test  is  increased.   It  is  shown  (Table  10)  that  when  d 
is  approximately  If,  N  is  8  at  P  =  .  95»  when  the  type  II 
error  Is  considered  the  required  W,  to  detect  d  =*  I4.  when 
P  =  .05,  $  =  .3,  is  19.   The  effects  of  power  are  very 
great. 

estimates  of  Beta  were  found  to  differ  from  week  to 
week.   Duncan's  multiple  range  test  shows  that  the  Beta 
estimate  found  on  July  2  was  significantly  higher  tnan  all 
others.   Estimates  collected  on  June  6,  July  9,  June  2^, 
and  July  16  did  not  differ. 

Homogeneity  of  tolerance  was  more  evident  on  July  2. 
Groups  tested  on  other  dates  showed  greater  heterogeneity 
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than  did  this  group.   Four  of  the  six  lots  tested  showed 
basically  the  same  degree  of  susceptibility  to  the  toxicant. 

SUMMARY 

Variation  in  response  of  a  laboratory-reared  KUN  strain 
of  insecticide  susceptible,  four-day-old,  male  houseflies 
to  malathion  applied  topically  was  studied. 

Several  statistical  methods  were  used  to  detect  and 
measure  variation  due  to  various  sources.   It  was  found 
that  there  were  two  sources  of  variation.   The  first  was 
heterogeneity  between  per  cent  mortalities.   This  was 
expected  and  its  significance  was  detected  by  Bartlett's 
test  of  homogeneity.   It  was  found  that  this  within  type 
of  variance  was  greatest  when  the  proportion  affected 
approached  0.  >  and  was  less  at  higher  proportions. 

Using  a  component  analysis  it  was  found  that  a  greater 
source  of  variation,  that  attributed  to  lot  to  lot  of 
houseflies,  was  also  present.   This  second  type  of  variation, 
among  lots,  must  be  considered  in  obtaining  reliable  esti- 
mates of  insecticidal  potency. 

The  effect  of  the  second  source  of  variation  upon 
sample  size  is  discussed.   A  randomized  complete  block 
design  is  suggested  to  determine  estimates  of  parameters. 
Days  (lots  of  insects)  represent  blocks  and  their  effects 
can  be  removed  to  gain  efficiency  over  a  completely  random 
design. 
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Estimates  of  the  I&dQ  and  Beta  were  determined  by  the 
probit  analysis.  Estimates  of  these  parameters  were  found 
to  differ  when  they  were  determined  from  different  lots  of 
houseflies.  The  pooled  estimate  of  the  i^^0  is  comparable 
to  estimates  found  by  other  workers  at  this  Station. 
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The  effects  of  variances  encountered  when  bloassays, 
involving  quantal  data,  are  used  to  determine  potency  of 
insecticides  were  studied. 

Variation  in  response  of  insects  to  insecticides  have 
been  reported  by  several  workers.   Davidow,  et  al.  ,  (1955) 
Yamasaki  and  Ishii,  (1955)  and  Sun  and  Sun  (1953)  described 
various  circumstances  known  to  cause  variation  in  response. 
Uunther,  et  al. ,  (1958)  found  greater  reproducibility 
within  repetitions  than  among  repetitions  done  on  different 
days;  they  used  Tribollum  confusum  Div.  as  an  indicator. 

The  KUN  strain,  male,  lj.-day-old  housefiies,  reared  for 
many  years  In  the  laboratory  according  to  the  standardized 
CSMA  procedure  and  treated  topically  with  malathion  carried 
in  acetone,  was  used  to  determine  response  to  the  insecti- 
cide, malathion.   ii.ach  experiment  was  conducted  each 
Tuesday  between  June  18  and  July  23,    1958. 

Several  statistical  methods  were  used  to  detect  and 
measure  variation  due  to  various  sources.   It  was  founu 
that  there  were  two  sources  of  variation.   The  first  was 
heterogeneity  between  per  cent  mortalities.   This  was 
expected  and  its  significance  was  detected  by  Bartlett's 
test  of  homogeneity.   It  was  found  that  this  within- type 
of  variance  was  greatest  when  the  proportion  affected 
approached  0.5  and  was  less  at  higher  proportions. 

Using  a  component  analysis  it  was  found  that  a  greater 
source  of  variation,  that  attributed  to  lot  to  lot  of 


houseflies,  was  also  present.   This  second  type  of  variation, 
among  lots,  must  be  considered  when  obtaining  reliable 
estimates  of  potency. 

The  effect  of  the  second  source  of  variation  upon 
sample  size  is  discussed.   A  randomized  complete  block 
design  is  suggested  to  determine  estimates  of  parameters. 
Days  (lots  of  insects)  represent  blocks  and  their  effects 
can  be  removed  to  gain  efficiency  over  a  completely  random 
design. 

Estimates  of  the  LIVq  and  Beta  were  determined  by  the 
probit  analysis.   Estimates  of  these  parameters  were  found 
to  differ  v.hen  they  were  determined  from  different  lets 
of  houseflies.   The  pooled  estimate  of  the  UV~  is  com- 
parable to  estimates  found  by  other  workers  at  this  Station. 


